Generalized Pearson distributions for charged particles 
interacting with an electric and/or a magnetic field 



A. Rossani a , A.M. Scarfone b 

a Istituto Nazionale per la Fisica della Materia ( CNISM-INFM) , Dipartimento di Fisica, 

Politecnico di Torino, 1-10129, Italy 
b Dipartimento di Fisica and Istituto Nazionale per la Fisica della Materia (CNR-INFM) , 
Sezione del Politecnico di Torino, 1-10129, Italy 



Abstract 

The linear Boltzmann equation for elastic and/or inelastic scattering is applied 
to derive the distribution function of a spatially homogeneous system of charged 
particles spreading in a host medium of two-level atoms and subjected to exter- 
nal electric and/or magnetic fields. We construct a Fokker-Planck approxima- 
tion to the kinetic equations and derive the most general class of distributions 
for the given problem by discussing in detail some physically meaningful cases. 
The equivalence with the transport theory of electrons in a phonon background 
is also discussed. 
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1. Introduction 

It is nowadays widely accepted that the Boltzmann equation constitutes 
a very powerful mathematical model to study different kinetic processes like 
fluid-dynamic problems, chemical and nuclear reactions, diffusion and others 
1] . Notwithstanding, the complicate mathematical structure embodied in this 
equation makes hard to obtain explicit solutions for specific problems, both from 
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an analytical or a numerical point of view. 

In Q], it has been developed a formalism to introduce consistently the inelas- 
tic interactions in the Boltzmann equation. In particular, it was considered a 
mixture of test particles (TP) which spread in a medium of field particles (FP) 
endowed with two levels of internal energy. In this way, under suitable hypothe- 
ses, a system describing the transport equation for the diffusion of TP in the 
medium has been derived. 

A particular simple assumption, but still preserving physical interest, that the 
mass M of FP is much greater than the mass m of TP, gives rise to a model in 
which the TP can gain or loose a fixed amount of energy. 

In this paper, we present a generalization of the results given in Q by consid- 
ering the case of charged TP diffusing in a medium (still made by FP endowed 
with two energy levels) , whose interactions arc modeled not only by means of an 
inelastic collision integral (as in Q) but also by the presence of elastic collision. 
In the most general fashion, we assume the presence of an external electric 
and/or magnetic field interacting with the system. 

In the case of TP interacting with FP only by means of elastic scattering, it is 
usual to adopt a Fokker-Planck approximation for the collision integral. Such 
approximation leads to a solvable equation for the distribution function (see 
[3, |j|). As it will be shown in this paper, the same approximation can be em- 
ployed also in presence of the inelastic scattering, which leads again to a solvable 
equation for the distribution function. 

Depending on the type of interactions occurring in the system, we obtain sev- 
eral physically relevant distributions. Among them, the generalized Pearson 

genau-Druyvesteyn distribution [4( are obtained 
in the case of hard sphere interactions, whilst power-law distributions Q, S] and 
modified power-law distributions |g| are derived in the case of Maxwellian in- 
teractions. 

Under the above mentioned hypotheses on the masses of TP and FP, the the- 
ory we are presenting can describe the diffusion of light particles in an heavy 
medium, so that both a loss and a gain of a fixed amount of energy is possible. 
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In this way, we can show how the present model is connected with the transport 
theory of electrons in a semiconductor lattice [llj, [llj . 

The paper is organized as follows. Section II is devoted to the description of 
the physical situation we deal with. In Section III, by employing the truncated 
spherical harmonic expansion (Pi approximation), we derive a system for the 
first two components N(v) and J(v) of the distribution function f(v). Then, 
by using the same method which is commonly applied to in the presence of 
elastic collisions only, we construct in Section IV a Fokker-Plank approximation 
for the general case containing both the elastic and the inelastic terms. Ex- 
plicit distributions corresponding to hard sphere and Maxwell interaction law, 
in the presence of electric and/or magnetic fields and elastic and/or inelastic 
scattering, are derived in Section V, whilst in Section VI, we obtain the expres- 
sion of the particles and heat currents related to some of these distributions. 
In Section VII, we study the proprieties of the inelastic collision integral only 
and discuss some aspects like the collision invariants and the trend to equilib- 
rium. Finally, in Section VIII, we show the mathematical equivalence of the 
Boltzmann equation without the elastic collisions and the transport theory of 
electrons in a phonon background. Conclusive comments are reported in Sec- 
tion IX, whilst an Appendix contains a detailed information on the distribution 
functions obtained in this paper. 

2. Outline of the problem 

Consider a spatially homogeneous medium of FP with mass M, endowed 
with one excited internal energy level. We call AE > the gap of internal 
energy between the excited and the fundamental level. Through this medium 
we consider TP, endowed with mass m and charge Q, which diffuse in the 
presence of an external electric field E and/or magnetic field B. The TP are 
supposed to interact with the medium according to the following scheme 

TP + FPi ^ TP + FP 2 , (2.1) 
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where FPi and FP 2 represent the fundamental and excited state of FP, whose 
number density will be denoted by N\ and A/2 , respectively. The number density 
n of TP is considered much lower than the number density Af = A/i + A/2 of FP, 
so that, the medium can be modeled as a fixed background in thcrmodynamical 
equilibrium at the temperature T. 
According to statistical mechanics we have 

N 2 ( AE\ 

= wn 



A^ exp b^J- (2 - 2) 

In this case, the kinetic equations for the distribution function of TP / = 
f(x, v, t) can be written as follows 



df df Q / \ df fdf\ cl (df 

dt dx m V ) dv \dtJ coU \dt /coU 



(2.3) 



The elastic collision integral is given by 

ol 



^ =J JgF\g,n-n')[f(v')F(w')-f(v)F(w)\dwdSl', (2.4) 



M \ 3/2 / Mu 1 



where 

T = M \j^f) ex H~i^J' (2 - 5) 

and I el (g, ft ■ ft') is the elastic cross section, g = \v — w\ is the relative speed 
with v and w the velocities of the incoming particles, whereas the post-collision 
velocities are given by 

v '= 2{m \ M ) (™v + Mw + Mgrt) , (2.6) 

w' = (mv + Mw-Mg ft') . (2.7) 

2(m + M) v y ' y ' 

Moreover, ft = (v — w)/g is the unit vector in the direction of the relative speed, 

whilst the two-dimensional unit sphere S 2 is the domain of integration for the 

unit vector 11'. 

Differently, for the inelastic collision integral we assume that M is much larger 
than m so that the velocity distribution function Tk{v) = ^fe(x, v, t) for FP^, 
at equilibrium, can be approximated by the still-particle distribution given by 

T k (v) =Af k 5(v) , (2.8) 
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that is a Maxwellian with M 



oo. 



In this way, the inelastic collision integr; 



al becomes 



( 



dj_ 

at 



) 



in 




jViv 2 + r n {v+, n-n')f(v+n') 



coll 




(2.9) 



where U(x) is the step function, v± — \Jv 2 ztrj with r\ = 2 AE/m and I ln (v, ft- 
ft') is the inelastic cross section of the TP-FPi collision. In Eq. (|2.9p . the first 
integral represents the gain of particles with velocity v due to the collisions of 
particles with velocity v± and unexcited/excited scatterers, respectively, whilst 
the second integral represents the loss of particles with velocity v due to colli- 
sions with unexcited/excited scatterers, respectively. 
In the following, we assume 



with a = 'el' or 'in', which allows to perform the Fokker-Planck approximation 
both in presence of the elastic and the inelastic collisions. These conditions are 
surely fulfilled for the hard sphere and the Maxwellian interactions examined in 
this work. 

3. Spherical harmonic expansion 

As usual ;Sf], if both the spatial gradients and the electric field are small 
enough, we may adopt the Pi approximation, which is, basically, a weighted 
residual method where the trial solution is given by a truncated spherical har- 
monic expansion for the distribution function of TP with shape functions 1 and 



\imv 2 i a {v, n-n') = , 



(2.10) 



ft: 



f(v) = N(v) + ft ■ J(v) , 



(3.1) 
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where 



and 



J(w) = — / . 
4tt J 



(3.2) 



(3.3) 



Starting from the Boltzmann equation (|2.3p . we derive the evolution equations 
for the unknowns iV(w) and J(v) by projecting it over the weights 1 and ft 
(Galerkin method) as follows 



£-N(v) + ~ V • J(v) + e ~ |- p J(«) 



97V 
"97 



coll 



97V 

ITT 



coll 

(3.4) 



p iv)+v VN(v) + eJL N( V ) + 6 x J(„) = ^ + ^ 

(3.5) 

where we set e = QE/m and 6 = QB/rn. 

The expression of the elastic collision integrals are well-known in literature and 
we refer to the relevant text-books (see, for instance, Ref. [3j). 
Instead, the inelastic collision integrals become 



ON 



dJ 

~dt 



1 r 



coll 



JV(w+)M v\ + N(v-)Af2 v 2 C(v) U{v 2 _) 



- N(v) [M 2 v\ I?(v + ) + M v 2 I?(v) U(v 2 _) 



1 r 



coll 



j(«+)M ^ ij>+) + j (v_ ) 7V 2 « 2 /j» c/(^) 

J(«) \N 2 v% /j>+) +M v 2 ff(y) U{v 2 _) 



where 



r e a (v) = 2tt / r n (v, fj,)fi l dfi , 



(3.6) 



(3.7) 



(3.8) 



with I = 0, 1 and fi = ft ■ $Y. 

In view of the manipulations we are going to perform, by assuming a stationary 
and space independent case, it is convenient to introduce the auxiliary quantities 
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£ = v 2 , F(g) = N(v) and G(£) = J(v), so that Eqs. (JUJ) and become 



(9 



<9F 



<9F 
^7 



3^ d£ 



coll V / coll 
el 



(3.9) 
(3.10) 



coll 



where 



OF 
~dt 



coll 



i r 

71 

+F(£-ri)Mtofr(£)U(£-Ti) 
F(0 



(3.11) 



dG 
~~dt 



coll 



l r 

71 

G 



(3.12) 



and erf (£) = 

4. The Fokker-Planck collision term 

In the following we extend a standard procedure, often employed to derive 
the Fokker-Planck approximation of [d f /dt) e to the case of inelastic colli- 
sions. 

We assume 

/hT\ 2 „ hT hT /hT\ 2 

(4-1) 



kT 

IT 



« 7] ~ — << 

M m 



m 



which means that rj is considered small with respect to the thermal square speed 
of a TP but still large with respect to the thermal square speed of a FP. 

Under this assumption we can adopt a procedure which leads to the Fokker- 
Planck approximation of (dF/dt)™^. In fact, by considering an arbitrary 
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smooth function $(£), one can easily show that 



We expand the integrand at the right hand side in power series of rj and retain 
only the terms up to r; 2 , so that 



oo 

o 

oo 


oc ✓ 



777 /'y 

2^(0 + ^(0 



where we integrated by parts. Since $(£) is arbitrary we can set 
'dF\ in N d 



777 /7 



(4.3) 



(4.4) 



and from the stationary condition (dF/dt) 1 ^ = we derive the following 
equilibrium distribution function for the present approximation 



F(0 = C exp 



2kT 

that is a Maxwellian with C the normalization constant. 

Let us now observe that Eq. (|3.12[) . together with Eq. (14. 4p . shows that G(£) 

C(t7 2 ) so that we can approximate (9G/di)™u in 



(4.5) 



coll 



where 



^(0 = ^(0-^(0 . 



(4.6) 



(4.7) 
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is the transport cross section for the inelastic process. 

Now, we want now utilize the approximate interaction terms (|4.4|) and (|4.6[) 
when both inelastic and elastic collisions occur. 

It is known [3] that for m « M the Fokker-Planck elastic collision terms reads 

cl 



OF 
~dt 



coll 



and 



_2_Af d_ 

M 7f af 



dG 



mF (0 + 2kT-^F(0 



where 



^ = -A^G(£K>(£) 

° 1 / coll 



(4.8) 

(4.9) 
(4.10) 



is the transport cross section for the elastic process. 

In particular, if cr ln (^) ~ c cl (0> ^ ^ s actually possible to add together the two 
collision terms. In fact, we can say that the order of approximation is the same 
if rj 2 ~ kT^/M ~ (k T) 2 / (M m) , which corresponds exactly to our assumption 
on r] 2 . 

Finally, by collecting the effects of elastic and inelastic interactions we obtain 
the following equations 



2 d_ 
3 6 ' dl 



^<(0 + ^ 2 - el (0' 



and 
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2 e ^ F(0 + -j= x G(0 = -G(0 <r{£) 



(4.11) 



(4.12) 



with <r(£) = <r in (£) + ct o1 (0, e = e/W and b = 6/A/\ 
We can solve this last equation for G(£) to obtain 

' e x b e b 



G(0 = -2 



<n0 + j 



6 + <r(£) e 



(4.13) 



that, introduced in Eq. (|4.11|) . gives 



*H0 + -f 



.M0 



a{0e 2 
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2kT 



no 



This equation is easily integrated, giving 



= C exp 



2fcT 



3(0 



where the constant C is obtained by the normalization condition 

oo 

2ir J y/iF({) d£ = n, 



(4.14) 
(4.15) 

(4.16) 
(4.17) 



and the function H(£) is given by 

s(0 = 



F(x) dx 



(4.18) 



with 



F(x) = I (^ 2 (x) + 6 2 ) 



a{x) 



+ e 2 x <j{x) 



+ 



AkT 



}f xo*(x) + -^o?(x) 



(4.19) 

The distribution (14.1611 contains a rich variety of cases, some of them known in 

nn 

the literature, like, for instance, the Margenau-Druyvesteyn distribution [3|, |4J 
the power-law distribution [3,3, the modified power-law distribution Q and, 
for e = 0, the Maxwell distribution. 

5. Some physically meaningful cases 

Under certain assumptions on the cross sections we can obtain the explicit 
expression of the distribution function (|4.16|) . In particular, we consider two 
physically meaningful interaction laws which give rise to different situations. 
The first case is given by the hard sphere interactions, by assuming constant 
values for all the cross sections, with er cl (£) = <r cl and cr in (£) = er m . Equation 
(|4.18p becomes 



m = J 



1 



a\ x + ci2 



(03 x + 04) (a 5 x + a 6 ) 



dx 



(5.1) 
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■p(0 = 0, (5.3) 
3spondi 



with ai = 8cre 2 /3, a 2 = 8 (e • b) 2 /(3cr), a 3 = 8kTa cl /M, a 4 = rf a„> a 5 = <? 2 
and ag = b 2 . 

This integral can be easily calculated, obtaining the distribution 

F 1 (0 = Ci (l + h £ + 6 2 C 2 ) exp ( - 6 4 € - 65 arctan(6 6 + 6 7 £)) , (5.2) 

which depends on the parameters b\, . . . , bj, whose relationships with the mi- 
croscopic quantities of the system are given in Appendix. 

Function (|5.2p is a generalized Pearson distribution, since it can be derived as 

a solution of the following linear differential equation 

dp(Q w 6 £ 2 + w 5 g + w 4 
d(, w 3 £ 2 + w 2 £ + wi 

for suitable values of parameters w\ , . . . , tog. The case z«6 = corresponds to 

the original problem defining the family of the Pearson distributions 

Several interesting distributions can be obtained, as special cases, from the 

family (|5.2|) . In particular, without the elastic collisions (a cl = and a = a" 1 ) 

the distribution reduces to 

F 2 (0 = C 2 (l + ci C) ^ exp (-c 3 , (5.4) 

where the coefficients Cj > (see Appendix) include the effects of the inelas- 
tic scattering only. This function is a modified power-law distribution. It has 
been previously considered in [9( from a kinetic approach, based on a non- 
linear Fokker-Planck equation taking into account the effects of a generalized 
exclusion-inclusion principle 

Another interesting case arises from Eq. (|5.2p without the magnetic field 

F 3 (0 = C 3 (l + d 1 ^j d2 exp (-rf 3 6 , (5-5) 

which is the Margenau-Druyvesteyn distribution with positive coefficients di > 
0, whose expression in terms of microscopical quantities are given in Appendix. 
We summarize in Table 1 the main distributions obtained from the hard sphere 
interaction law in the presence of elastic and/or inelastic collisions and in the 
presence of electric and/or magnetic fields. 
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Table 1. Classification of the distribution functions for the hard sphere 
interaction case with elastic and/or inelastic collisions and electric and/or 

magnetic external fields. 





a cl and a in 




er in 


e and b 


Generalized-Pearson 


Generalized-Pearson 


Modified Tsallis 


e 


Druyvesteyn 


Druyvesteyn 


Heated Maxwellian 


b 


Maxwellian 


Maxwellian 


Maxwellian 



The second case we study is given by the Maxwellian interactions where cross 
sections are given by c el (£) = cr el /y/£ and cr m (£) = ^"VV?- Formula (|4.18[) now 
becomes 

fix 



m = 1 1 1 





dx 



f2X + f 3 

with ft = 8 [(e • bf/a + ae 2 ]/[i (a 2 + b 2 )}, f 2 = 8 kT5 cX /M and f 3 = rj 2 < 
By performing the integration we obtain 



(5.6) 



F 5 (£)=C 5 1+91 a exp(- ff3 



-92 



(5.7) 

which is again a modified power-law distribution, depending on three positive 
coefficients gi > (see Appendix). This distribution is functionally equivalent to 
the distribution (|5.4[) , i.e. hard sphere TP with purely inelastically interacting 
behaviors like Maxwellian TP. 

An interesting situation, already derived in [l^ , is obtained when the scattering 
is purely inelastic, corresponding to g$ — > limit. This gives rise to a power-law 
distribution that can be written in the form 

I 1/(1-9) 

2kT\ 

known in literature as Tsallis distribution [l^ . 



F e (0 = C 6 



!-(!-<?) 



(5.8) 



As show in Appendix, the de- 
formed parameter q is related to the microscopic quantities of the system ac- 
cording to 

M\ hT (f ■ h\ 2 4- p 2 n 2 

(5.9) 



^_ i+ 16/cT (e-b) 2 + e 2 a 2 



3 to (b 2 -\- d' 2 ) Tj 2 o (Tq 11 
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which is an increasing function of the electric field and a decreasing function 
of the magnetic field. (Remark that in Ref. 12 1 a different definition for the 
deformed parameter q has been adopted, which is related to the present one 
through q — > 2 — q) . 

In Table 2 we summarize the main distributions obtained from the Maxwell 
interaction law in presence of elastic and/or inelastic collisions and in presence 
of electric and/or magnetic fields. 

Table 2. Classification of the distribution functions in the Maxwell interaction 
case with elastic and/or inelastic collisions and electric and/or magnetic 

external fields. 





a cl and a in 


-el 


-in 


e and b 


Modified Tsallis 


Heated Maxwellian 


Tsallis 


e 


Modified Tsallis 


Heated Maxwellian 


Tsallis 


b 


Maxwellian 


Maxwellian 


Maxwellian 



6. Physical quantities 

From a physical point of view, we are interested to study the conditions 
assuring the existence of the first £ momenta of the distribution. They are 
defined in 

, fv'f(v)dv 
{V) Jf(v)dv ' (bA) 
so that, for £ = 1, we obtain the definition of the density current, given by 



(6.2) 



(6.3) 



with u the mean velocity of TP while, for £ = 2 we have 

(v ) = - r , 

which corresponds to the mean-squares speed, where, for sake of convenience, 
we introduced the quantity 



c t = J rVKK. 




(6.4) 
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with 

#0=eKp(-£Ls(O) . (6.5) 

The existence of these quantities is assured for almost all the distributions de- 
rived in the previous section, since the presence of the exponential factor guar- 
antees a fast convergence of the integrals. 

An exception is given by the distribution (|5.8j) . which has a purely power-law 
asymptotic behavior. In this case, the convergence of the above integrals, and 
in particular the normalization condition, imposes a limitation on the value that 
the parameter q can assumes. We recall, (see [15| and reference therein), that 
in the stationary and homogenous case the transport equation for charged par- 
ticles subjected to an external electric field has not, in general, a normalizable 
solution. The nonexistence of such a solution is known as runaway phenomenon 
and its occurrence depends on how fast is the decay of the cross section as a 
function of the energy. 

For the distribution (|5.8p , the convergence of the quantity (v 2 ) (which also 
guarantees the convergence of the lower momenta with £ < 2) is assured by the 
further condition 

1 < q < \ , (6.6) 
5 

which defines, through Eq. (|5.9p . a boundary in e and b for the onset of the 
runaway phenomenon. 

Another physically relevant quantity is the heat flux, defined by 
q = -to J (v — u) (v — it) 2 f(v) dv 

= n m\Q - jj ^ + u 2 \ u , (6.7) 



where we have set 



oo 

Q = ^U 2 G{i)di. (6.8) 
6 n u J 



In general, both the quantities j and q cannot be obtained in a closed analytical 
form. 

Notwithstanding, exact expressions for these currents can be obtained for the 
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case of charged TP governed by the only inelastic scattering. A simpler assump- 
tion, still preserving the physical interest, is given by setting 6 = 0. In this case, 
our problem shrunk in just one dimension, since now both the currents j and 
q assume the same direction of the electric field. In this simple situation, for 
particles undergoing to hard sphere interactions, we obtain the expressions 



Q = 



3 n 
n m 



2 7r k T* 
2 k T% 



1/2 



1/2 



2 m e 



3 y/rr \ m J 7r V 3 fc T* a 
where T*, derived in Appendix, is given by 

T. = T ( 1 



- -1 



3 r] 2 a o-J, 11 



(6.9) 
(6.10) 

(6.11) 



Let us observe that in the e — > limit the particle current shows an Ohmic 
behavior with j cx e, while for e — > oo it approaches an asymptotic value 

In the same way, in the case of Maxwell interactions, we have 



e 

3 =n- 



and 



q = nm — e , 



(6.13) 



a <j J 
and the particle current shows always an Ohmic behavior. 
More in general, it is known that in the limit of weak electric field, particle 
current and heat current can be expressed by means of the following constitutive 
equations 



J 

q 

where, the Einstein law 



= kuE + fci2 VT + k 13 Vn 
= kn E + k 22 VT + k 23 Vn 



fcn 



k 2 3 
k 2 i 



can be rewritten in the compact form 
j = kn E + k 12 VT , and 



kT 

n\Q\ ' 



q = k 21 E + k 22 VT 



(6.14) 
(6.15) 

(6.16) 

(6.17) 
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where E = E — {kT/n\Q\)V n. 

Further, we can extract informations about all the coefficients fcy from the Eqs. 
(|6.2p and (|6.7p although these relations hold only for the spatially homogeneous 
case. 

In fact, as known, the coefficients kij are related to each other by means of the 
following relations: 



k\2 

that is the Onsager law and [16| 



-T , (6.18) 



that is the Weidemann-Franz law 
Starting from the expression of ku and &2i derived from the linearization of j 
and q, we can obtain all the other coefficients. For instance, we can linearize 
the expression of the currents (|6.9|) and (|6.10|) , which hold for the hard sphere 
interactions, to obtain the coefficients 

11 = o — k? r' rm " ) ' anc i ^21 



3 ctTV ' ; ' 3 o-AT 



7r m 



(6.20) 

from which one can deduce the expressions of the remaining coefficients k^ . 



7. Study of the inelastic collision integral 

In this section, we deal with the properties and the consequences that the 
inelastic collision integral induces into the system. For this purpose, we start 
by considering the weak form of the transport equation, without the elastic 
collisions. This can be accomplished by considering a smooth function <f>(v) and 
by introducing the following auxiliary functional 

0M - /*»(§{)> 

n')-M /(»«)" 

xv 2 + I in (v + , n ■ n')vdv . (7.1) 



dfl / dfl' 



(vsi) - 4>{v+ n') 
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From this equation, we find that any arbitrary function 4>{v) — <£>(u 2 ) (depend- 
ing only on the modulo of v) and periodic with period e 2 



$(v 2 + e 2 ) = $(w 2 ) , 



(7.2) 



implies G[<f>] = 0, i.e. it is a collisional invariant for the given problem. 
In particular, the choice $ = 1, a constant, corresponds to the TP number 
conservation. We observe that the initial assumption M — > oo has the important 
consequence that TP does not conserve the ordinary physical quantities like 
energy and momentum. Moreover, Eq. (|7.2p implies the existence of infinite 
constants of motion (spurious invariants), where $ = 1 is just the simplest case. 
Let us now consider the above weak equation applied to the function 

„2 1 



(v) = In 



/(«) exp 



mv 
~kT 



(7.3) 



In this case, the auxiliary function G[(f>] becomes 

g[<t>] = JdflJdSl'J {in /(««)- In [£/(«+«')]} 

£ f(v+ ft') - f(v ft)] JV 2 v\ I in {v+, n-n')vdv < . (7.4) 

Equation (|7.4| is a non positive quantity, as it follows from the inequality (a — 
b) ln(a/6) > and consequently it states an H theorem for the present problem. 
At equilibrium, it must be Q[4>] = 0, for any v, Q and Q' . This implies the 
following relation 

\nf(vn)-]n[£f(v + n')]=0, (7.5) 



from which one easily obtains a family of equilibrium distributions 

,,, i' 2 

f(v) = T(v 2 ) exp 



2kT 



where T(v 2 ) is an arbitrary function that fulfils the periodicity condition 

T(v 2 + e 2 ) =T(v 2 ) , 



(7.6) 



(7.7) 



originated from the condition l|7.2[) . 

Clearly, the family of distributions (14. 16)) contains, as a special case, for r(i> 2 ) a 
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constant, the Maxwellian distribution, whilst, in all the other cases, it generates 
a family of infinite distorted Maxwellians, one for each choice of the function 



8. Equivalence with transport of electrons in a phonon medium 

In order to better clarify the role of the Boltzmann equation with the inelastic 
collisional integral only, we shall show now that this equation is mathematically 
equivalent to the transport equation for electrons in a phonon background. In- 
deed, the results obtained in the previous section are surprisingly similar with 
those found by Majorana in the field of electrons transport in a semiconductor 




ialui- 



First of all, let us rewrite Eq. ([13]) with (df/dt) 



cl 



0, in the following form 



CO. 



,11 




(8.1) 




QB/(me) and f(w) 




+Af 2 w 2 r n (w, si ■ si') /(«;_ n') u(wi) dsi' 




+M w 2 i in (w, si • si') u{ w 2 _) dsi' , 



(8.2) 



where w± = Vw 2 ± 1 and I in (w, SI ■ Si') = el in {v, SI -Si'). 
By taking into account the relation 



U(w ± 1) S(w' - w±) =2w± 8(W 




(8.3) 



we may express (df/dt) 



in 



in the following equivalent form 



coll 



18 



£K(w, w', n ■ n') 5(w' 2 - w 2 - 1) 



where 



+fC(w', id, Ct ■ n') (5(id 2 - id' - 1) /(«;') du/ 
-/(w) y fC(w, id', f2 ■ n') <j(u/ 2 - ui 2 - 1) 

+ejc(v/, w, n • n') <$(u/ 2 - w 2 - i)l <w , 



/co, u/, n ■ n') = 27V 2 — r n (w', n ■ n') 



(8.4) 



(8.5) 



In Ref. jlol |. free electrons have been considered interacting with monochro- 
matic phonons (energy tko) of a semiconductor lattice. The kinetic equation for 
electrons reported in that paper turns out to have exactly the same form as Eq. 
(|8.1|) . when the last expression for (d f /dt)™ oll is adopted, where AE plays the 
role of Hlo. 



9. Conclusions 

Starting from the Boltzmann picture, we have studied a system of charged 
particles spreading in a medium made by two-level atoms in the presence of 
an external electric and/or magnetic field. For the two meaningful interaction 
laws of hard sphere and Maxwellian particles we derived two wide families of 
distribution functions. 

Concerning the hard sphere interactions, while the appearance of generalized 
Pearson distribution is fascinating, we must summarize the assumptions we 
utilized to derive this result: 1) the assumption that M is much greater than 
77i, which leads to the linear equation for inelastic scattering; 2) the assumption 
that the electric field is small, so that the Pi approximation is allowed; 3) the 
assumption that AE is much smaller than the thermal energy of a TP, so that 
the Fokker-Planck approximation is allowed; 4) the assumption that the cross 
sections are governed by power-laws. 

On the other hand, concerning the Maxwellian interactions, we would like to 
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remark that the statistical parameter q, for the power-law distribution in the 
Tsallis form, can actually be derived strictly from microscopic arguments. This 
circumstance, encountered in few models in which the Tsallis distribution arises, 
is a very important point, related to the physical validity of this distribution, 
which is often argument of intensive debate [ItI [3| . 

APPENDIX. 

In this appendix we give some details about the distributions derived in the 
section V. 



Hard sphere interaction 

The hard sphere interaction follow by assuming constant values for all cross 
sections, with ct o1 (£) = <r cl and cr in (£) = er in . In this case, Eq. flUTS]) takes the 
form of a generalized Pearson distribution 

Fi(0 = Ci (l + h £ + b 2 C 2 ) cxp ( - 64 £ - 65 arctan(6 6 + b 7 £)) , (A.l) 



whose coefficients b\ , 
system according to 

61 = 

62 = 

63 = 

64 = 

b 5 ■■ 

be ■■ 

b 7 - 

with 



67, are related to the microscopical parameters of the 



8 
12 



kT a l & 



2 „el 



M (3 2 ' 
{kT) 2 <ra cl 
mM e 2 



2fcT 

m 
2kT 
Pi 



M e 2 /3i 16 (e • bf 
3kT aa el ~ 1 a 



1 



16 



kT a 2 a el 



M /3 3 ' 



kT 



(A.2) 
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8 (e ■ b) 2 2 2 in 

P2 = - h V CT 

O (T 



03 - 



32^a 2 a cl / 3 2 -/? 2 



1/2 



The following particular cases arise: 

1) For cr cl = 0, we obtain the modified power-law distribution 



^2(0=^2(1 + 01^) ^ cxp(-c 3 £) 



with 



ci = 



C2 



C3 = 



6 2 



1 + 



8 (e-bf 



3 rj 2 b 2 a a$ 



4 to e 2 6 2 - (e ■ b) 



3 fcT 7 7 77 2 cr 3 ^ n 

TO 



where 



7?\ 



and 



7 = 1 + 



8e 2 



2) For b — 0, we obtain the Margenau-Druyvesteyn distribution 



F 3 {0 = C 3 (l + d 1 t) 2 exp(-rf 3 , 



with 



M e 2 + l^craj, 11 ' 



d 2 
d 3 



6(fcT) 2 crcr cl ' 
m 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



2kT ' 

3) Again, for b = but without the elastic collisions (cr cl = 0), we obtain the 
heated-Maxwellian distribution 

to£ 



with 



F 4 (0 = C 4 exp 



T* = T 1 



2 & 



8e 2 



3 7] (7 CTq" 



(A.9) 



(A.10) 
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Maxwellian interaction 

The Maxwellian interaction follow by assuming <r el (£) = cr cl /^/t; and cr m (£) = 
er 111 / v^- The function (|4.16|) takes the form of a modified power-law distribution 

F 5 (0^C 5 (l+ gi ^y 92 exp(-<feO , (A.11) 

whose coefficients <7i, (72 and 33 are related to the microscopical parameters of 
the system in 

„ k T s a cl 



,9i 



If (Jq 



" = A"****- < A - i2 » 



.93 



2kT s 



where 



M (e ■ b) 2 + e 2 a 2 



The following particular cases arise: 

1) For (j cl = and a — a m , we obtain the power-law distribution 

F 6 (0 = C 6 (l + h^y h2 , (A.14) 

with 

8 (e ■ b) 2 +e 2 5 2 , , to . 
«i = ^— — - — — , and ho = . A. 15 

1 3 n 2 (b 2 + a 2 )aa% 1 2kTh 1 V 1 

Distribution (|ATl4|) can be written in the form (|5.8|) . by introducing the de- 
formed parameter q = 1 + I//12, given by 

3 to rj z (b z + a z ) a 

2) For ct 111 = and <r = cr , we obtain the heated Maxwellian distribution 



where 



r f = (i + Or, and C^ y^/ . (A.18) 
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